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We describe a novel type of magnetic domain wall which, in contrast to Bloch or Neel walls, is 
non-localized and, in a certain temperature range, non-monotonic. The wall appears as the mean- 
field solution of the two-dimensional ferromagnetic Ising model frustrated by the long-ranged dipolar 
interaction. We provide experimental evidence of this wall delocalization in the stripe-domain phase 
of perpendicularly magnetized ultrathin magnetic films. In agreement with experimental results, 
we find that the stripe width decreases with increasing temperature and approaches a finite value 
at the Curie-temperature following a power law. The same kind of wall and a similar temperature 
dependence of the stripe width is expected in the mean-field approximation of the two-dimensional 
Coulomb frustrated Ising ferromagnet. 
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Introduction. — Recent rigorous works [l|, 0] state that 
the spontaneous magnetization of a two-dimensional (2d) 
Ising ferromagnet vanishes exactly if a dipolar coupling 
is present (Dipolar Frustrated Ising Ferromagnet, DFIF). 
Frustrating interactions on different spatial scales occur 
e.g. in ultrathin magnetic films where the spins point per- 
pendicular to the film plane P, i, i, i, i, 0, MM, 0, [HI, 
[l2 . 13, 14|. Several approaches indicate a striped ground 
state [Jll, H, [H, [111 ■^itl^ spi^ modulation along one 
in-plane direction and uniform spin alignment along the 
orthogonal in-plane direction. At finite temperatures, 
such spin "microemulsions" llTII suffer from the Landau- 
Peierls instability [J, i, i, EllE [13, III [11 that delocal- 
izes, in the thermodynamic limit, the stripe positionjthus 
reducing the positional order to quasi-long range @,[^,[^. 
Notice, however, that the stripe width remains well de- 
fined at finite temperatures even in the thermodynamic 
limit [l8| , and that experiments on real systems do indeed 
observe the persistence of stripe order at finite tempera- 
tures [1,0, 11, 12, [ill, possibly because of domain- wall 
pinning [3, [20, [21 1 . It is thus worthwhile to study the 
Mean-Field (MF) behavior of a DFIF at finite tempera- 
ture in the hope that some characteristics are "robust" 
enough to be valid beyond the MF approximation and 
to be observable in experiments. Competing interactions 
acting on different length scales are fundamental to many 
different chemical and physical systems [1, [13, [H, [H, [25| 
so that robust MF results on such a general model like 
the DFIF may have a wide significance. 

A central question that motivated this work is the equi- 
librium stripe width h* (number of lattice parameters) 
at finite temperatures. One result appears to be well 
established in 2d: The stripe width in the ground state 
depends exponentially on the ratio between the exchange 
(J) and the dipolar (g) energy per atom [H. [sl. [ol. [l6l. [l7| . 
The temperature dependence, instead, is controversial. 
Within the MF approximation, h* is expected to decrease 
with temperature, because it should reach a finite value 



on the order of ^ at the temperature Tc of the MF tran- 
sition to the paramagnetic state where the spin averages 
to zero at every site [S 
on sharp interfaces [26 



15 1 . Theoretical arguments based 
23, [Hi predict a (stretched) expo- 



nential decrease of h* (T) down to an atomic length scale 
at the transition temperature Tc. Within the spherical 
approximation, the modulation length of a related model 
(the Coulomb Frustrated Ferromagnet) "monotonically 
increases with temperature until it diverges at a disorder 
line temperature" , and this independently of the dimen- 
sionality of the system [2^ . Experimental results on the 
temperature dependence of the stripe width are contro- 
versial as well [13, [13, [13 and in some cases experiments 
do not show any change of h* with temperature [^ . 

In spite of its apparent simplicity, a detailed study of 
the stripe width of a DFIF as a function of temperature 
within the MF approximation has not been reported yet. 
Here, wc fill this gap and solve the relevant MF equa- 
tions for the DFIF model on a discrete lattice, finding 
a number of unexpected results. 1) The sharp-interface 
assumption gives the correct stripe width at low temper- 
atures, but fails to reproduce the results of the full MF 
calculation close to the transition temperature Tc. 2) 
Near Tc, the temperature dependence of h* crosses over 
to a power law. 3) This cross-over is accompanied by a 
delocalization of the wall between adjacent stripes: The 
profile changes from square-like at low temperatures to 
cosine-like at Tc. 4) At intermediate temperatures, the 
interface between two adjacent stripes develops a pro- 
nounced non-monotonic "shoulder" tailing down toward 
the center of the stripe according to a power law. The 
non-local walls are in striking contrast to the profile of 
the domain walls encountered in the Ising model without 
dipolar interaction [s^l and also with the shape of con- 
ventional Bloch or Neel walls dividing domains in typical 
Heisenberg or planar ferromagnets 31[ . The strength of 
the "shoulder" structure at intermediate temperatures 
depends on the relative strength of J and g, but it oc- 
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curs over mesoscopic scales and in a sizeable range of 
temperatures so that it should be observable by spatially 
resolved experiments which have a high enough signal- 
to-noise ratio. The cosine-like profile, instead, is realized 
sufhciently close to Tc independently of the ratio j. Here 
we provide experimental evidence that the spin profile of 
the stripes changes indeed from square-like at low tem- 
peratures to cosine-like at T^. 

The model. — The DFIF Hamiltonian on a discrete lat- 



tice reads TL — 



N 



where 



i and j are two-dimensional indices (i = {ix^iy}) with 
1 < ia; < -^a: and \ < iy < Ly, (7i — ±1, means that 
the sum is restricted to nearest neighbors, {i ^ j} indi- 
cates a sum over all the different pairs in the lattice, and 
N — Lx ■ Ly is the total number of spins. This system 
is treated in the MF approximation and the spin super- 
structure is assumed to be a stripe pattern (a very com- 
mon pattern in real systems), i.e. we require the averages 
"^(ix.jy) ^ {'^(i^,iy)) translationally invariant along 

the y-direction and periodic along the x-direction with a 
period of 2h: rn(^i^^i^) = rm^ and rrii^+ah = 
where a G Z. The relevant MF quantities are thus re- 
duced to the h independent variables rm^. To deter- 
mine the equilibrium configurations, we first find the pro- 
file for given h from the self-consistent MF equations 



= tanh(/3E,.V; 



where /3 



V/^ is the effective interaction matrix [35 
imize the free energy per spin 



l/fceT and 
then we min- 



f 



2h ^ ^ 



h 



with respect to h. 

The stripe width. — Typically, in thin magnetic films 
J ^ g. We are thus interested in large ratios S = J/g, 
for which we recover the value of h*{T — 0) known from 
Ref. [11] and the finite value h*{Tc) known from Ref. 0. 
Tc and the limiting magnetization profile rui^ [T — > Tc) 
(see Fig. 2a) are found as the maximum eigenvalue and 
the corresponding eigenvector of V^^ j ^ . 

In Fig. 1, we plot h*{T) for 5 = 10. The polygonal 
appearance of the graphs is due to h* changing only by 
±1 in a discrete model. In order to understand the low- 
temperature behavior of h*{T) and the low-T plateau, 
we introduce a manageable "sharp-interface" two-spin 
model, where mi = rrih = rriw and m^^ = for 
ix (lower left inset of Fig. 1). In the limit of 

large J, rrib is just the MF value for a ferromagnetic 
Ising model, m?, = tanh(4/3 Jm?,), and is the corre- 
sponding MF value for a spin adjacent to a domain wall, 
rriyj = tanh(/3 J(mu, +mb)). Due to the reduced exchange 
energy in the argument of the tanh, mw{T) < mb{T) 
(see later Fig. 2a). Inserting mw{T) and mb{T) into the 
sharp-interface free energy and minimizing it with re- 
spect to /i, we obtain /i2spin('^) 



exp(l-|-A/4(7m^), where 
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FIG. 1: Equilibrium stripe width. The solid line represents 
h* (T) as a function of the reduced temperature T/Tc for 5 = 
10. The dots on the line are those points for which we show 
spin profiles in Fig. 2a. The upper inset plots h* (T) — h* (Tc) 
versus {Tc - T)/Tc in the critical region for 5 = 100. Full 
(empty) dots correspond to the upper (lower) corners of the 
steps of the /i*(T)-curve (see sketch, lower right). A fit to 
the full dots is shown. The fit of the full (empty) dots gives 
an exponent of 2.0 (1.9). The dashed line shows ^2spin(T) 
calculated within the two-spin model (sharp-interface limit) 
for S = 10. The lower left inset is a cartoon visualizing the 
two-spin model. 
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The dashed curve in Fig. 1 representing /i2spin(^) fo'^ 
(5 = 10 reproduces the low-temperature behavior of 
the numerical solution but fails at higher temperatures, 
where it gives ^2spin(^c) ~ 1. In the upper inset of Fig. 1, 
we plot h*{T) - h* {Tj) vs (Tc - T)/Tc for the fuU MF 
calculation close to Tc in a log-log plot, showing that 
the domain width behaves asymptotically according to a 
power law h*{T) ~ h*{Tc) ~ (Tc - T)^, with A « 2 as 
discussed in the figure caption. This numerical outcome 
seems to confirm the conjecture of Ref. |19|, but is at 
odds with the sharp-interface limit of Refs. [26l. [27l [28| . 
It would be interesting to review the experimental results 
of Ref. [H, [2^ under the point of view of a power law. 
The argument of Ref. [l^ associates the cross-over to the 
power-law behavior with the higher harmonics (responsi- 
ble for the sharp interface at low temperatures) vanishing 
with increasing temperature and thus with the broaden- 
ing of the spin profile to a cosine-like profile close to Tc. 

The magnetization profile. — The spin profiles m^^ at 
different temperatures, obtained from the transcenden- 
tal MF equations, are plotted in Fig. 2a for selected 
temperatures (marked with dots in Fig. 1). We iden- 
tify three regimes: (i) A low-T regime, corresponding 
to the plateau in the h*{T) curves, with a square- like 
profile, (ii) An intermediate T regime, corresponding to 
the steep descent of h*{T). Here, a novel feature con- 
sisting of a double-shoulder and a wall delocalization are 
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FIG. 2: (a) Spin profile within a stripe calculated for S — 
10 at different temperatures and the corresponding values of 
h*{T). The triangles locate the interface spins, (b) Spin 
profile for a single domain wall with 5 = 10 (solid line) and 
5 = oo ("Landau" profile, dashed line), for T /Tc = 0.95. The 
magnetization approaches the m(cxD)-value exponentially for 
(5 = 00, but it decays as 1/ix for finite 5 (lower inset). The 
upper inset shows the demagnetizing field (dots) originating 
from the dipolar interaction for a step-like spin profile (solid 
line). 



observed, (iii) A high-T regime, corresponding to the 
critical region, where the magnetization has indeed the 
cosine-like profile also expected from analytical consider- 
ations and leading to the power-law behavior of the equi- 
librium stripe width [T9| . Notice that only within the 
first regime is the interface sharp and does the bulk mag- 
netization 7712, • ■ • J f^h-i (circles) not change very much, 
so that the two-spin model (dashed curve in Fig. 1) is 
indeed justified. In order to understand the origin of 
the non-monotonic shoulder and the wall delocalization 
we have solved the MF equations for a single domain 
wall. In Fig. 2b we compare the profiles in the absence 
(dashed line) and in the presence (solid line) of the dipo- 
lar interaction. For 5 = oo {g — 0), the profile is of the 
Landau-type; it increases monotonically and attains the 
asymptotic value exponentially. For finite 5, the profile is 
not monotonic: it has a shoulder close to the wall center 
and then decays to m(-l-oo) as the inverse of the distance 
(lower inset). Both features derive from the dipolar (de- 





FIG. 3: (a) SEMPA [11, [SJ measurement of a stripe sec- 
tion at low temperature (10 K). The spin polarization is en- 
coded by a gray scale. The measured image is 4000 pixels 
(17 /^m) wide and 5 pixels (21 nm) high. For better inspec- 
tion, the image has been stretched by a factor of 70 along the 
vertical direction. Within the spatial resolution of the exper- 
iment, this mesurement corresponds to 5 scans of the same 
scan line, displayed as successive lines in one image. Ther- 
mal drift causes some displacement between the scan lines, 
(b) SEMPA measurement of a stripe section at high temper- 
ature (330 K). The measured image is 400 pixels (4.25 /im) 
wide and 35 pixels (370 nm) high. This image is displayed 
in its real proportions, (c) Experimental spin profiles across 
one stripe at low temperature (10 K, h=9 fim, empty circles) 
and close to the transition temperature (330 K, h = 430 nm, 
black squares) as extracted from images (a) and (b) respec- 
tively. To obtain the profile, the scan lines are aligned to 
compensate for the thermal drift and then the aligned scan 
lines are averaged. To further reduce the noise level, these 
raw profiles are resampled by averaging a number of adja- 
cent points. Both profiles have been normalized to the same 
width and amplitude. The horizontal error bars represent the 
spatial resolution (±50 nm) of the experiment. It is deter- 
mined from the topographic profile of a sharp edge, obtained 
by the same procedure from a topographic image measured 
with exactly the same parameters as the magnetic image 3b. 
For the low-T profile, the horizontal error is on the order of 
the size of the data points. The vertical error bars are due 
to the statistical noise of the secondary electrons counted in 
SEMPA [igl . [3^ . Note that no traces of an in-plane compo- 
nent of the spin polarization can be found in simultaneously 
recorded images of the in-plane spin polarization. Details of 
the experiment will be published in a more extended review 
and are available on request. 



magnetizing) field i7'*'P(ia;), which is plotted for a step- 
like profile in the upper inset. Close to the center of 
the wall, H'^^P{ix) almost vanishes because of the com- 
pensation between the fields generated by up and down 
spins. There, the deviation from the Landau-type wall 
(dashed line) is small. The approach to the _ff'^'P(-|-oo) 
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value occurs as 1 ji^ and depresses the spin profile below 
the asymptotic value for g = 0. We have also checked 
within a continuum model that the demagnetizing field 
far from the wall vanishes in the infinite-thickness limit 
so that the three-dimensional, monotonic and localized 
"Landau" wall is recovered in 3d. Thus, the formation 
of the non-monotonic long-ranged wall is a purely two- 
dimensional effect. 

The square-like profile at low temperatures delocalizes 
into a cosine-like profile close to Tcj no matter how small 
the dipolar interaction is. We provide experimental 
evidence of this wall delocalization in Fig. 3. The spin 
profile was measured in SEMPA ^] experiments on 
ultrathin Fe films grown epitaxially on Cu(lOO). These 
systems are magnetized perpendicularly to the film plane 
and show the sought-for stripe structure [IJ. The two 
different profiles at low temperature (empty dots) and 
close to the stripe-paramagnetic transition temperature 
(full dots) point to the realization of the MF cross-over 
shown in Fig. 2a. 

Conclusions. — In summary, we have shown that when 
dealing with "spin microemulsions" (and probably with 
analogous pattern- forming systems), the range of valid- 
ity of the sharp-interface limit must be evaluated care- 
fully and that important physical features, like the stripe 



width, crucially depend on whether the actual interface 
is sharp or not. In addition, we have discovered that the 
"Landau" -type walls (and probably the Bloch- and Neel- 
type walls as well) must be modified in low-dimensional 
systems because the dipolar interaction produces a non- 
monotonic, long-range tail which is absent in 3d systems 
such as those considered by Landau [lH. Our study 
has focused on the DFIF model on a discrete lattice. 
However, its results appear to be relevant [111 for the 
Coulomb Frustrated Ising Ferromagnet [s^ as well. In 
the Coulomb system, antiferromagnetic interactions de- 
cay as rather than |r as the dipolar interaction, 
but preliminary results indicate that the phenomenology 
(domain shrinking, power-law approach to a finite value 
at Tc, delocalization and presence of a shoulder in the 
domain wall profile) is the same. Although our work 
is based on the MF approximation, it produces results 
which appear to be realized in real pattern-forming sys- 
tems, such as the power-law dependence of the stripe 
width [l^ and the spin profile (Fig. 3c), and it might 
provide a reasonable starting point for more sophisticated 
theoretical work. 
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